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We study the effect of a half wave rectified sinusoidal electromagnetic (EM) wave on the Kitaev 
honeycomb model with an additional magneto-electric coupling term arising due to induced polar¬ 
ization of the bonds. Within the framework of Floquet analysis, we show that merging of a pair of 
Dirac points in the gapless region of the Kitaev model leading to a semi-Dirac spectrum is indeed 
possible by externally varying the amplitude and the phase of the EM field. 


I. INTRODUCTION 

The Kitaev honeycomb model [T] is an exactly solv¬ 
able anisotropic spin-^ model with Ising-like interactions 
cr^cr^/, and assigned to the three bonds in 

the honeycomb lattice; stand for the a{= x, y, z)- 

th Pauli spin operator residing on the r/r'-th site of the 
lattice. Although there has not been an experimental re¬ 
alization of the Kitaev honeycomb spin model till date, 
there exist some promising proposals of experimental re¬ 
alization of the spin models involving ultracold atoms [2] 
and polar molecules [3] in optical lattices. There have 
also been theoretical proposals that the Kitaev model 
may be realized as an effective low-energy model for a 
Mott insulator with strong spin-orbit coupling [UlS]. The 
model exhibits a rich phase diagram where depending on 
the relative strength of coupling between the Ising inter¬ 
actions the system either has a graphene [6] like gapless 
spectrum or a gapped spectrum; the gapless phase and 
the gapped phases are separated by the phase bound¬ 
aries where the spectrum is of semi-Dirac nature laiHi 
such that the associated quantum critical point [nmni is 
of an anisotropic nature m- 

In the gapless phase of the model the dipersion relation 
e{k) vanishes at the contact points between two bands. 
At these so called Dirac points the spectrum is linear and 
due to time reversal symmetry of the Hamiltonian these 
points occur in pairs in the reciprocal space at points D 
and —D. It has been observed that for graphene-like sys¬ 
tems varying the hopping parameters, by application of 
periodic perturbations, leads to a shift in the point loca¬ 
tion and may lead to merging of two such Dirac points. 
On merger the resulting spectrum is a semi-Dirac spec¬ 
trum with a dispersion relation which is quadratic along 
one direction in the reciprocal space and linear along the 
other (perpendicular) direction. The merging of Dirac 
points is accompanied by a topological phase transition 
from a semimetallic to an insulating phase [IMS]. We 
note that the magnetic field dependence of Landau lev¬ 
els in a Graphene-like structure resulting in a semi-Dirac 
spectrum was already reported much earlier [nHin]. 

It is worthwhile to mention here that the study of topo¬ 
logical insulators and the topological protection of edge 
states [2Qll45] have gained a huge importance in recent 


years. In parallel, an extensive amount of work have 
also been dedicated to study in depth, driven closed 
quantum systems from the perspective of dynamical 
saturation [46]. dynamical localization [T7ll5Q] . dynamical 
freezing [51]. dynamical fidelity [52]. defect pro duct ion [53l 
l54] and thermalization[55|. But it is the advent of irradi¬ 
ated (Floquet) grapehene[44l l56l |57j and Floquet topo¬ 
logical insulators that have gelled the two, facilitating 
experimentalists to be able to probe and verify a few of 
the many recent theoretical work done on the effect of 
driving on closed topological quantum systems |40H42]. 

Recently, a large body of work has also been done 
in which a time-periodic perturbation has been used 
to dynamically change the phases of Hamiltonians 
exhibiting different topological properties [HI |20l [43]- 
|45l (56] [58H6T] . In the context of Graphene, Delplace 
et al US] used an external time-periodic electromagnetic 
(EM) field to merge a pair of Dirac points in the high 
frequency regime by varying the amplitude and the 
phase of the field applied. On the other hand, in our 
work we look at the equivalent problem in the Kitaev 
model. However, there is a subtle issue that needs to 
be addressed; the Kitaev model unlike graphene consists 
of localized spins which do not possess linear momenta 
and hence the minimal coupling scheme used in m 
cannot be extended to our case. To overcome this 
situation, we employ the scheme introduced by Sato 
et al [62] where the effect of an external EM field on the 
Kitaev honeycomb model was studied by incorporating a 
magneto-electric coupling term in the Hamiltonian. We 
use this modified Hamiltonian to study the possibility 
of merging the Dirac points in the gapless phase of the 
original Kitaev model through an external eliptically 
polarized, half wave rectified electromagnetic wave. Such 
a manipulation allows us to control the phases and hence 
the topological properties of the Kitaev model driven by 
an external periodic perturbation in the high-frequency 
limit. 

The main motivation of our work is to externally tune 
the couplings of the Kitaev Hamiltonian in such a way so 
that the Dirac points in the gapless region merge lead¬ 
ing to a semi-Dirac spectrum within the gapless region 
of the undriven Kitaev model itself. In our subsequent 
discussion we will present the fact that unlike Graphene 
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with nearest neighbor hopping, the linear dispersion (or 
the Dirac point) in the gapless region of the undriven Ki- 
taev model is anisotropic (with different velocities in the 
two different directions). Moreover, the presence of lo¬ 
calized spins on lattice sites brings forth further difficulty 
while trying to couple the spins to any perturbation asso¬ 
ciated with an external electric field, in stark contrast to 
the case of Graphene, where the electron momentum can 
easily couple with the applied electromagnetic field via 
Peierls’ substitution. Hence, in spite of all such difficult 
issues how merging of Dirac points in the gapless region 
of the undriven Kitaev model can still be achieved via 
a periodic external electromagnetic field, is an intriguing 
question of non-equilibrium statistical mechanics. The 
fact that the Kitaev model is the only exactly solvable 
model in two dimensions provides us with further moti¬ 
vation. 

The rest of the paper is organized as follows; the Ki¬ 
taev model is outlined in SecHIl We describe the mod¬ 
ifications necessary to make the Hamiltonian couple to 
an EM wave in Sec |HI[ T his is followed by a review of 
Floquet theory in Sec |IV| generalizing it to the context of 
the space periodic Kitaev model in SecjV| We discuss the 
effect of driving on the Hamiltonian in Sec |VI| and the re¬ 
sulting phase diagram in SecjVHj Concluding comments 
are presented in Sec |VHl| 


II. KITAEV MODEL 


The Kitaev model consists of spin-l/2’s placed on the 
sites of a honeycomb lattice with a Hamiltonian of the 
form: 

-f^Kitaev = + 

j-\-l=even 

+ ( 1 ) 

where j, I are the column and row indices respectively, 
^ are Pauli matrices, representing the spins, at the 
site labeled (m,n) and J^’s are the corresponding cou¬ 
pling parameters (see Fig. Q. In this section the energy 
spectrum for the time independent Kitaev model and the 
corresponding phase will be reviewed before we proceed 
to discuss how the addition of an EM coupling between 
the spins and an external field changes the spectrum by 
affecting the interaction strengths (Ji’s). We set Ji > 0 
below without any loss of generality. 

Referring to Fig{^ we take the unit cells of the lat¬ 
tice to be the vertical bonds with sites labeled A and 
B, respectively so that there are N/2 unit cells for a 
system of N sites. Choosing the nearest neighbor dis¬ 
tance to be l/\/3 we label each unit cell by a vector 

n = irii + ^ 2 : with ni ^2 denoting the coor¬ 

dinates of the B site in the unit cell n. The spanning 
vectors joining the neighbouring unit cells are then given 



FIG. 1: (color online) This shows the honeycomb lattice of 
the Kitaev model with the interactions labeled as Ji. The vec¬ 
tors Ml ,2 span the lattice, and the positions of the sublattice 
points are labeled by A and B. 


by Ml = -h and M 2 = ~ To fermionize 

the system we define the Majorana operators as: 



for j 1 = even 
for j 1 = even 
for j = odd 
for j 1 = odd 


( 2 ) 


These are Hermitian operators representing real fermions 
(i.e., = a) satisfying the anti-commutation relations 

for example, CLm',n} = In terms of 

these operators the Hamiltonian Q takes the form: 

^ ~ ^ “h T ‘ (3) 

ft 


Remarkably, operators Dfi = b'^a'^A defined on each pla- 
quette commute with each other and with the Hamil¬ 
tonian and their eigenvalues can take the values ±1 in¬ 
dependently for each n This enables us to reduce 2^- 
dimensional Hilbert space into 2^/^ dimensions. How¬ 
ever, it has been established that the global ground state 
of the model lies in the sector in which Dfi = 1 for all ft. 
To extract the spectrum we move into the Fourier space 
summing over half the Brillouin zone due to the Majo¬ 
rana nature of the fermions. The Hamiltonian can then 
be written in the (2 x 2) Bogoliubov-de Gennes form as: 


Hk= 2 


Ji sin(fc.Mi) — J 2 sin(fc.M 2 ) 


( 4 ) 


+2 


J 3 + Ji cos(/c.Mi) + J 2 cos(/c.M 2 ) 
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FIG. 2: (color online) This shows the phase diagram of the 
Kitaev model. The plane Ji + J2 + J3 = 1 in the parameter 
space described by (Ji, J2, ^/s) has been chosen for illustration 
without any loss of generality. The phases marked ^1,2,3 are 
gapped while the phase marked B is gapless. We highlight 
the two phases (A and B), that we have concentrated on, in 
the discussion. 


where are the Pauli matrices denoting the pseu¬ 
dospins. The dispersion relation can then be immediately 
derived as 

= 4 sin(^.Mi) — J 2 sin(^.M 2 )^ (5) 

+4 + Ji cos(£Mi) + J 2 cos(fc.M2)^ . 

The phase diagram of the model as deduced from Eq. (§ 
is shown in Fig. Given that > 0, it is convenient to 
choose a normalization Ji = 1 which describes points 
lying within (or on) an equilateral triangle. This triangle 
can be divided into four smaller equilateral triangles. 
They are labeled as Ai where Ji > J 2 + Js^ ^2 where 
J 2 > Ji + J 3 and A 3 where J 3 > Ji + J 2 - In the region 
B each of the Ji is less than the sum of the other two. 
It turns out the spectrum is gapped in all the A regions, 
and gapless in the B region . 

We reiterate that the purpose of our work is to tune 
the couplings of the Kitaev Hamiltonian given in Eq. 0 
through an external control, such that Dirac points in 
the B region merge leading to a semi-Dirac spectrum 
even within the gapless phase of the unperturbed Kitaev 
model. This we achieve using a periodic external elec¬ 
tromagnetic field where the phase and amplitude of the 
field renormalizes the hopping strengths J^’s leading to 
the desired merging. 


III. LIGHT IMPINGEMENT ON THE KITAEV 
MODEL 

As discussed above that the Kitaev model is essentially 
a system of localized spins on a lattice without any 
real momentum which could couple to the EM field as 
happens in the case of itinerant electrons in Graphene. 


To surmount this difficulty, Sato et al [62] introduced a 
magneto-electric coupling between the EM field and the 
spin lattice. 

In our approach, it is assumed that the EM field affects 
the underlying lattice by causing the stretching (or com¬ 
pression) of the underlying lattice bonds those host the 
spins; therefore, the EM fields leads to an effective polar¬ 
ization enabling us to define an effective polarization ten¬ 
sor that couples to the external field; this coupling conse¬ 
quently alters the bare bond strengths {Ji, J 2 ^ J 3 ) renor¬ 
malizing them to {J'l^ J 21 Foi* clarity, we would like to 
stress that the co-ordinate axes 1, 2 and 3, are oriented in 
the bond-direction in the spin space while (x, y, z) rep¬ 
resent spatial cartesian coordinates. The ME coupling 
originates from electric polarization on each bond and its 
strength is assumed to be proportional to the exchange 
interaction of the bond [62|. The Hamiltonian now in¬ 
cludes an additional piece 


H{t) = i^Kitaev + ^ME ( 6 ) 


One assumes a polarization of the form P(r,r')c — 
Hc^cr^cr^,, defined on each bond (r, r')^, where Hq, is the 
ME coupling vector and a = 1,2,3 denote the direction 
of the bonds. This leads to a total polarization as Aotai = 

Ea E(r,r')a HME{t) = -E{t) ■ Ptotal- We 

have taken the polarization along the bonds of the Ki¬ 
taev model which lie in a plane of the Honeycomb lattice 
which in the Cartesian basis {x^y) is represented as: 


n = n3y+till (x - Vsy) + Ilia(-:r - Vsy). (7) 

In the subsequent discussion, we shall work with a 
Hamiltonian of the form: 


H = 


E 


/' T'^l _i_ T '^2 2 I t /^3 ^3 \ 


j-\-l=even 

(8) 

using an electric field of the form E{t) = 
A(sin(f2t), sin(Dt + 0),O)), where the renormalized 
coupling strengths are = {J^ — H^.E) and the index 
a runs over the bonds (1, 2, 3). 

Our work does rely on the magneto-electric (ME) cou¬ 
pling term (first introduced in Ref. [62] in the case of 
the Kitaev model) to produce a merging of the Dirac 
points by light impingement; we would, however, like 
to emphasize that such an interaction term is not un¬ 
common in multiferroic materials (which exhibit both 
ferroelectric and magnetic ordering) and has been ex¬ 
perimentally detected in rare earth metal (R = Tb, Gd) 
perovskite manganites [RMn 03 ] [63]. RMn 03 are hexag¬ 
onal polar crystals with ferrimagnetic, ferromagnetic or 
anti-ferromagnetic orders [64]) and may be implemented 
in optical lattices as proposed in Ref. [3] to generate the 
Kitaev spin model. The ME coupling is strong in such 
multiferroic materials due to the presence of frustrated 
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magnetic exchange interactions and non-collinear spin or¬ 
der. Theoretically, the ME coupling term can easily be 
understood (in terms of the phenomenological theory de¬ 
veloped by Landau) by expanding the free energy as a 
function of the magnetization and the electric polariza¬ 
tion. As a result the linear ME effect arises from the 
cross term where P^y is a second ranked ten¬ 

sor, and Py are the magnetization and the electric 
polarization vectors, respectively. This coupling also de¬ 
pends on the symmetry of the crystal and is only present 
when both the time reversal symmetry and the inver¬ 
sion symmetry are broken. We recall that the breaking 
of time reversal symmetry and the inversion symmetry 
are necessary for spontaneous magnetization and ferro¬ 
electric polarization, respectively. We refer to Tokura et 
al. [65] for a discussion of the exchange-strict ion type 
of ME coupling used in our paper (in particular in the 
discussion below Eq. [6]) and of the ME coupling in mul- 
tiferroics. In our model, we have assumed that such a 
coupling originates from a slight distortion of the under¬ 
lying lattice due to the electric field of the laser light or 
through a tiny phonon-mediated dimerization leading to 
a polarization of the bonds whose strength as a result is 
some fraction of the strength of spin-spin exchange inter¬ 
action along the bonds. 

IV. FLOQUET-BLOCH THEORY IN 

PERIODICALLY DRIVEN SYSTEMS: 

The Eloquet technique [66ll68] (which is a tempo¬ 
ral version of Bloch’s theorem) deals with Hamiltoni¬ 
ans subjected to a time-periodic potential of the form 
H = iLo + E(t), where V{t P r) = V{t). Recalling the 
discrete time translation operator T defined through the 
relation + r) = Aq,' 0 q,(x, t), we note 

that for stationary solution, Aq, has to be a pure phase of 
the form . Thus we have a solution of the form : 

'lpa{x,t + T) = (9) 

where Ua{x, t P r) is a time-periodic function that satis¬ 
fies the condition Ua{x,t P r) = Ua{x^t) and uja = 

In a spirit similar to Bloch theorem, where one defines 
quasi momenta, in the time periodic case one introduces 
the notion of quasi energies of the form 
defined within the first Brillouin zone, (“^5 ^)- 
When viewed stroboscopically (at the end of each 
complete period r), the problem effectively reduces to a 
time-independent problem; the unitary time evolution 
operator after n-complete periods assumes a simple form 
given by I/(nr, 0 ) = T exp (—i = [t/(r, O)]’^, 

where T denotes the time ordering operator. Working 
in a representation in which U is diagonal such that 
Ud = diagonal[e~'^^^], one can evaluate the quasi 
energies in a straightforward manner. 

In the present case we are dealing with a Hamiltonian 
that is periodic in both space and time and is represented 


by: 

H{x P di,t P r) = H{x,t), (10) 

where a^’s are the lattice vectors. The Eloquet-Bloch 
Hamiltonian is then given as: 

•H(r, k) = e-**-*iJ(r)e**-* - idr 

n{T, k) = H^{t) - idr (11) 

where is in the A?-space. The Eloquet states, de¬ 

fined earlier, now get generalized to Eloquet-Bloch states 
which are periodic in both space and time and obey the 
equation (i?g(r) - idr) g); where g) 

is periodic in both space and time. The ket is defined 
in a composed Hilbert space (the Sambe space) which is 
the direct product space of the original Hilbert space and 
and the space of time(r)-periodic functions [69]. The in¬ 
ner product in this space is defined by a composed scalar 
product: 

<< ... >>= - [ < ...> dt (12) 

^ Jo 

where < ... > is the normal inner product in Hilbert 
space. 

The above discussion leads to the conclusion that us¬ 
ing the Eloquet-Bloch ansatz, one can reduce both k 
and time to parameters in a time independent eigenvalue 
equation although the dimensionality of the Hilbert space 
gets augmented. This reduction enables us to exploit the 
nature of the quasi-energy spectrum to study the prop¬ 
erties of a driven Kitaev chain. Here, we have multiple 
copies of the same Bloch bands in frequency space, known 
as Eloquet Bloch bands, while the coupling between the 
copies is determined by the frequency of driving. We 
shall work in the high frequency limit where there is no 
coupling between the different frequency Brillouin zones 
restricting the system to the original Bloch band with 
renormalized hopping parameters. 

V. FLOQUET-BLOCH HAMILTONIAN FOR 
THE DRIVEN KITAEV MODEL 

In the case of the driven system, given in Eq. (§ , the 
couplings are time dependent. Consequently, to diago¬ 
nalize this system we will need to introduce time depen¬ 
dent Majorana operators (a, b) which lead to the form 

H = i jt) + 

ft 

+J'^Dfi(t)bfi{t)afi{t). (13) 

Using the time dependent form of Eq. @ in Majorana 
representation with the notation a ^ Ca and 6 —> C 5 , we 
obtain a time dependent (2 x 2) Hamiltonian of the form: 

a, (3 


( 14 ) 
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with both a, [3 running over the sublattice index (a, 6). 
The Hamiltonian matrix is given as: 


H = 


0 p{k,t) 
p*{k,t) 0 


(15) 


with p{k,t) = and 

Ar =2 J[{t) sin(^.Mi) — J^it) sm{k.M 2 ) 


(16) 


=2 

k 


J^{t) + J[{t) COs(^.Mi) + J^it) COs(^.M2) 


As the system is time-periodic the Majorana operators 
can be decomposed into their frequency Fourier modes 
through g(t) = The Hamiltonian 

given in Eq.(|14|) goes to: 

"jW = E E ixr- (17) 

p,p' a,/3 

Substituting the above relation in the Hamiltonian in 
Eq. 0 and recalling the definition in Eq.([T^, we ob¬ 
tain the following inner product, 


(('^a,k,p'\'^k(^)\'^l3,k,p)) ~ ^p'^p P^^p,p'^cx,/3 ( 18 ) 


with 


= 

- 


^ (19) 


Eq.(18) reveals that the static tight binding (TB) model 
in N dimensions under the application of an AC elec¬ 
tric field gets mapped to a time-independent TB model 
in {N + 1) dimensions with an effective static DC elec¬ 
tric field, i.e., the driving frequency (cj), acting along 
the added dimension. This effective electric field breaks 
the translational symmetry along the additional dimen¬ 
sion; we can now separate out two regimes on the ba¬ 
sis of high or low effective electric field, i.e, in terms of 
high or low driving frequency. Of course, one needs to 
compare it against the other energy scale in our model 
which is the tunnelling strength or the strength of the 
exchange interactions (Ji, 2 , 3 )- The N dimensional lat¬ 
tice thus repeats itself along the (asymmetrical with re¬ 
spect to frequency) extra dimension hosting lattice points 
that are linked with their neighboring lattice points via 
renormalised (or photon dressed) tunneling or exchange 
interaction strengths. 


VI. MERGING OF DIRAC POINTS 

To drive the system we choose a half wave rectified 
elliptically polarised electric field given by: 


E{t) 

E{t) 


T 

= A(sin(co’t), sin(co’t -h 0), 0), 0 < t < — (20) 


= 0 


T 

~2 


<t<T 



FIG. 3 : The Brillouin Zone (BZ) of the honeycomb lattice is 
shown, indicating the four TRIM points as has been stated in 
the text. 


where A is the magnitude of the field, the phase (j) is the 
measure of the ellipticity and the half wave rectified elec¬ 
tric field that repeats itself over a time period T with the 
frequency denoted by u; = The motivation behind 
using a half wave pulse (and not a complete sinusoid) is 
to prevent the integral in Eq. (18) from generating inter¬ 
actions of 5 function form rendering the dynamics trivial. 

Although the reason behind choosing a half wave rec¬ 
tified sinusoidal electromagnetic wave is to have an ef¬ 
fective renormalised Hamiltonian with competing terms 
of finite strength, these waves are indeed feasible exper¬ 
imentally and are commonly generated in optical setups 
via a Mach Zender optical modulator. This modulator 
functions by splitting the incoming electromagnetic beam 
into two waveguides where a time varying voltage in one 
of the waveguides induces a time-dependent phase shift 
(via the electro-optic effect) in the passing wave. When 
the undisturbed wave passing through one waveguide is 
made to interfere with the phase shifted wave from the 
other waveguide, the two waves interfere and the time- 
dependent phase shift manifests itself as an (time vary¬ 
ing) amplitude modulation of the initial wave. Thus, 
the amplitude of a laser beam may be tuned to obtain 
a half wave rectified sinusoidal electromagnetic wave by 
suitably varying the applied voltage (periodically) in the 
setup discussed above. 

One should now work in the high frequency regime 
(cj ^ Ti, J 2 , Js) to show that the above mentioned driv¬ 
ing scheme can make the Dirac points merge. The high 
frequency regime as has been previously mentioned, has 
been identified by comparing the two energy scales in 
the problem : the driving frequency and the strength of 
the exchange interactions. Now Eq. (18) tells us that 
the Eloquet Hamiltonian matrix is infinite dimensional 
in the space defined by as they belong to the set 
of all integers. Thus, at such high frequencies, the last 
term in Eq. (18) is more dominant, and hence the Elo¬ 
quet Hamiltonian matrix is roughly block diagonal. The 
Eloquet sub-bands are hence, uncoupled which allows us 
to set p ^ p' in Eq. (18). Hence we can restrict ourselves 
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to a system described by a (2 x 2) Hamiltonian similar 
to Eq.(p!^. This leads to a renormalization of the hop- 
J[ (Eq.p^). Eor our analysis we 


ping strengths Ji 
restrict ourselves to the block with p = 0. This implies 
that the quasi-energy eQ, = \p\. The system is studied 
using the spectrum generated for this reduced Hamilto¬ 
nian with modified coupling strengths recalling Eq. 0 
and the form of the electric field ( [2Q| ) (see Appendix for 
details) 


— J\ T — ^—a/S “h cos Hi 
A. / \ 

J 2 = J 2 ~\ -+ cos (pj YI 2 

2A 

J 3 = J 3 -Hs cosp (21) 

TT 

It is now necessary to obtain a condition involving the 
renormalized parameter , J 2 and J 3 for the merging 
of Dirac points. Due to time reversal symmetry of the 
Hamiltonian, the structure of the Dirac points in the k 
space is such that if there occurs a Dirac point at D, 
another Dirac point must be present at — D. When two 
such Dirac points merge, it implies that a point in the 
/c-space is mapped to itself under time reversal symme¬ 
try, this is called a time reversal invariant momentum 
(TRIM) point. To obtain the TRIMs for the Kitaev 
model we note that the two Dirac cones must move equal 
distances in /c-space along the lattice vectors Mi and 
M 2 . Thus giving the TRIM points as ^(pMi + qM 2 ) 
with (p, q) = ( 1 , 1 ), ( 1 , 0 ), ( 0 , 1 ), ( 0 , 0 ) (see Fig. @). This 
yields the following conditions on the hopping parame¬ 
ters: 

Mq —y J2 Tt/3 = 0 M\ —y = J2 

M2 ^ = + J3 M3 ^ 4 = 4 + 4 (22) 

These relations are the same as the ones for phase 
transition derived for the Kitaev model without the 
ME coupling (Eq.Q) [T]. The noticeable difference is 
now the conditions are imposed on the renormalized 
parameters J's. It should also be noted that with the 
renormalized coupling parameters, it is now possible to 
merge the Dirac points at the center of the BZ even if 
Ji ^ 0] this is because the merging conditions are 
imposed on the renormalized strengths J'’s and does no 
longer on bare coupling strengths J^s. The merging of 
Dirac points at the TRIM points produces a semi-Dirac 
kind of a spectrum which means that the spectrum is 
linear along one direction while being quadratic along 
the other. We can now vary A and p to achieve a 
transition from a Dirac to a semi-Dirac spectrum. 


VII. THE PHASE DIAGRAMS 

The phase diagram for the merging scenario will 
depend on all the parameters J^’s, A, p and the H^’s 



<f> 


FIG. 4 : (color online) The phase diagram for the anisotropic 
case 71 7^ 72 7^ 73- The diagram shows the lines in 

( 0 , A) space where the Dirac points merge. The position of 
merge is indicated by the Mi labels (Mo = blue (dotted) , 
Ml = green (solid) , M2 = brown (dot-dashed), M3 = pink 
(dashed)). The points where two such merge lines intersect 
represent regions where one of the couplings J- is zero. The 
phases are distinguished by the colour of the boundary walls 
as elaborated in the text. 



D 1 2 3 4 5 « 




FIG. 5 : (color online) The phase diagram for the isotropic 
case 71 =72 =73- The diagram shows the lines in ( 0 , A) 
space where the Dirac points merge. The position of merge is 
indicated by the Mi labels (Mi = brown (solid), M2 = pink 
(dotted), M3 = blue (dashed)). The condition for merge at 
Mo is missing because of the isotropy. The phases are distin¬ 
guished by the colour of the boundary walls as elaborated in 
the text. 
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(c) Semi-Dirac spectrum : The PDFs have now merged at one (d) Gapped quasi-energy spectrum at A = 2.1 

TRIM point (M2 here) at the toroidal BZ boundary at 
A = 1.382528373335 

FIG. 6: This k-space quasi-energy spectrum (denoted by E here) shows the dynamic merging of the PDFs in the gapless phase 
of the undriven Kitaev model for system parameters Ji — J 2 — J 2 , — \ with a field amplitude varying from (a) 0 to (b) 1.1 
going up to (c) 1.382528373335 and then to (d) 2.1 at a fixed phase of 0 = 1 and all 7^ = 0.5 (isotropic case). We see in (a) 
that there are 2 conical (anisotropic linear dispersion) PDFs touching the zero quasi-energy plane within the gapless phase of 
the undriven Kitaev model in the first BZ. As the field strength varies, we notice that although the quasi-energy spectrum is 
still gapless in (b), the 2 PDFs have shifted towards the BZ boundary. In (c), the amplitude of the field is such that the 2 
PDFs now lie at the BZ boundary. But, since the quasi-momentum (kx,ky) lie on a torus, the two PDFs are actually at the 
same quasi-momentum (TRIM) point. Hence, the Dirac points have merged to produce a semi-Dirac spectrum in the gapless 
phase of the undriven Kitaev model under driving at this TRIM point. We further observe in (d) that a slight variation of the 
field amplitude only manages to gap out the spectrum as expected. 


where only A and (j) serve as our externally tunable. We 
will next demonstrate how a pair of Dirac points (PDFs) 
located at the center of the undriven gapless region 
in the parameter space of J^’s can be merged through 
driving to produce a semi-Dirac spectrum. Two kinds of 
phase diagrams emerge depending on the anisotropies in 
polarisation. We will gradually discuss them below. 

To draw the phase diagrams we assume that the EM 
field affects the underlying lattice in such a way so that 
the polarizability constants (H^’s) will only be a fraction 
of the exchange constants (Ji’s). Thus, we can write : 

Ili = 7 i Ji (23) 

where 7 ^ < 1 . 


Let us now discuss the first scenario in which all 
the 7 i’s are unequal. The Fig.Q shows the different 
phases that can emerge from such a condition. The 
Dirac points may merge at a particular TRIM point 
Mi for different values of amplitude A and (j). Thus, 
the four continuous lines in four different colors (one 
color for each i = 0,1,2,3) in Fig.(|^ are the 
merging lines where a PDP get annihilated rendering a 
semi-Dirac spectrum. There are four such lines because 
there are four TRIM points as has been discussed earlier. 

We can go from one phase to the other in the phase 
diagram only upon crossing the merging lines. The 
phases, as a result, can be identified by observing 
whether in that phase a PDP can merge at all the four 
TRIM {Mi) points or it can merge at all the TRIM 


































points (M^’s) except one Mi where it can never merge. 
The phases in which the four different merging lines 
serve as the phase boundary are phases in which a PDF 
can merge at any one of the TRIM (Mi) points. There 
are also phases in which only three (or two) different 
merging lines form the phase boundary. In such phases, 
a PDP can can never merge at those Mi points (one or 
two respectively) whose merging lines do not constitute 
the phase boundary. The phases have thus, accordingly 
been identified in the figure (Fig. @). 

A secondary route towards merging is suggested by 
the structure of the renomalized J’s (see Eq.([^), 
namely the isotropic case where 71 =72 = 73 - The 
resulting phase diagram is shown in Fig.([^. The point 
to be noted in such a case is that the PDFs can merge 
at two or three TRIM points but not at all four. A 
pair of PDF’s can never merge at the central TRIM 
point (Mo) given by the condition + J 2 + J 3 = 0 , 
as is evident from the fact that for the isotropic case 
Si J'i — Si 7^ 0. 

This fact is reflected in the phase diagram. Character¬ 
izing the phases using the same classification scheme as 
used for the earlier anisotropic case; we find the existence 
of phases for which the boundary walls represent merger 
of two or three different TRIM points. It can also be 
seen that the phase represented by low values of the field 
amplitude A has three different boundary walls so from 
this phase we can effectively merge the PDF’s at all the 
TRIM points except at Mq. 

The phase diagram of the isotropic case (71 =72 = 73 ) 
is indeed very simple to understand. A plot of the 
quasi-energy spectrum {E) against the quasi-moment urns 
given by (kx^ky) as a means to illustrate the isotropic 
case is shown in figure (Fig. We vary the field 

amplitude A from 0 to 2.1 keeping the phase fixed at 
0 = 1. The quasi-energy spectrum changes from having 
a PDFs to a merger at a TRIM point on the M 2 line 
(pink line in Fig. (§) yielding a semi-Dirac spectrum at 
A = 1.382528373335. 

The Dirac points (DPs) which are time-reversed part¬ 
ner of each other can merge at the time reversal invariant 
momentum (TRIM) points only. There are seven TRIM 
points with one of them at the centre and the other six at 
the zone boundaries of a hexagonal Brillouin zone (BZ). 
But, since the TRIM points at the zone boundaries are 
equivalent modulo (or related via) a reciprocal lattice 
vector, there are only three inequivalent TRIM points at 
the zone edges of a hexagonal BZ. This takes the total 
tally of inequivalent TRIM points within a hexagonal BZ 
to four (three at the BZ boundary and one at the centre). 
To make the scenario further transparent, we refer to the 
Fig. 0 that explicitly tracks how the DPs undergo a 
merging transition as the vector potential (A) is varied. 
In the figure, we depict how a merging transition hap¬ 
pens when one DP moves and settles at one of the TRIM 
points while its time reversed partner settles at another 



b 

FIG. 7: (color online) The BZ boundary has been plotted 
in blue while the TRIM points lying on it are shown in hol¬ 
low red dots. We mark the DPs as little solid dots of vary¬ 
ing colors. There are only two DPs for a fixed vector po¬ 
tential strength (A) and hence, we see pairs of DPs in the 
same colour (and size) on the two opposite edges of the BZ. 
The growing size of the dots indicate the direction in which 
the DPs proceed nearly along the BZ boundary towards the 
TRIM points before merging when A has been varied from 
0 (DPs in small blue at the BZ corner which lies on the kx 
axis) to 0.2 (DPs in big red dot) to 0.4 (DPs in bigger cyan) 
to 0.6 (DPs in even bigger green) to 0.8 (DPs in large or¬ 
ange) to 1.0 (DPs in large magenta) to 1.2 (DPs in larger 
brown) to finally at A = 1.382528373335. The DPs merge for 
A = 1.382528373335 as they lie on top of the red TRIM point 
(black solid huge dot inside). The black arrow depicts that 
the two TRIM points in which the DPs merge differ only via 
a reciprocal lattice vector. 


TRIM point, where both the TRIM points are related by 
a reciprocal lattice vector. Here, one must note that the 
merging of DPs does not happen at the centre of the BZ; 
rather it occurs at two points at opposite edges related 
via a reciprocal lattice vector and hence are equivalent. 

A similar plot for the anisotropic case (71 7 ^ 72 7 ^ 73 ) 
can also be easily constructed. A careful analysis then 
would clearly corroborate the fact that the anisotropic 
case only differs from the isotropic case in the presence 
of the possibility that the PDFs can also come together 
at the Mo TRIM points. 


VIII. CONCLUSION 

In summary, through this work it is clear that exter¬ 
nal controls can be used to effectively move a PDP in 
the reciprocal space for the Kitaev honeycomb Hamilto¬ 
nian. In our analysis we have used a modified form of 
the Hamiltonian with an effective magneto-electric cou¬ 
pling to facilitate our manipulation of the location of the 
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PDF’s through an externally applied half wave rectified 
electric field. The effective phase diagram of the system is 
then critically dependent on the amplitude A and phase 
(j) of the external field. 

The important point is that through the scheme outlined 
the phases of the model can be externally controlled and 
indeed fine tuned phases can be created where the model 
reduces to an effective Ising model (at the intersection of 
two merge lines). Finally, we note that merging of DPs 
involves transition from gappless to gapped phases which 
may be experimentally captured in the entanglement en¬ 
tropy CO] as also theoretically predicted in m- 


of the time dependence in the Hamiltonian using the pe¬ 
riodic nature of the driving by diagonalising the Hamilto¬ 
nian in the Sambe space using the composed inner prod¬ 
uct as given in Eq. (p!8|).The further use of Eqs. ([T^ and 
(19 to perform the composed inner product yields : 


K'i = (Mt) - iBuit)) (A.l) 


where each of the can be cast as a 2 x 2 matrix of 
the form : 
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Hp'—p — 


0 Pp! —p{h^ 

Pp-pik) 0 


and pp'-pik) = 


(A.2) 


Appendix: The renormalised coupling constants 

In this appendix, we briefly discuss how the renor- 
malised coupling constants given by J's in Eq. (21) can 


be obtained. The main aim of this calculation is to get rid 


Since our interest lies in the high frequency regime 
where the Eloquet sub bands decouple and form 2x2 
blocks for each p, the limit p' ^ p is taken in the above 
Eq. (A.l) alongwith the choice of the p = 0 block as has 
been mentioned earlier, to obtain po(^) = ^o(^) “'^^ 0(^)5 
where ^o(^) and Hq(^) are : 


A',{k) 

B'oik) 


= 2 

= 2 


I Ji + — \/3 + cos0^ Hi I sin(^.Mi) — | J 2 + — + cos0^ n 2 | sin(^.M 2 ) 

{ 


(A.3) 


2 A 

J 3 - 113 cos 

TT 


i>| -h I Ji + — \/3 + cos0^ Hi I cos(fc.Mi) + I J 2 + — + cos0) n 2 \ cos(/c.M 2 ) 


Thus, it can easily be seen that in the high frequency 
regime, when our decoupled Eloquet sub blocks have 
the same form as our original Kitaev Hamiltonian (see 


Eq. (p!^), the terms within the curly braces are nothing 
but the renormalised coupling constants J's as has been 
provided in Eq. (21). 
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